The HLA region in ANCA-associated vasculitis: characterisation of genetic associations in a Scandinavian patient population.

OBJECTIVE: The antineutrophil cytoplasmic antibody (ANCA)-associated vasculitides (AAV) are inflammatory disorders with ANCA autoantibodies recognising either proteinase 3 (PR3-AAV) or myeloperoxidase (MPO-AAV). PR3-AAV and MPO-AAV have been associated with distinct loci in the human leucocyte antigen (HLA) region. While the association between MPO-AAV and HLA has been well characterised in East Asian populations where MPO-AAV is more common, studies in populations of European descent are limited. The aim of this study was to thoroughly characterise associations to the HLA region in Scandinavian patients with PR3-AAV as well as MPO-AAV. METHODS: Genotypes of single-nucleotide polymorphisms (SNPs) located in the HLA region were extracted from a targeted exome-sequencing dataset comprising Scandinavian AAV cases and controls. Classical HLA alleles were called using xHLA. After quality control, association analyses were performed of a joint SNP/classical HLA allele dataset for cases with PR3-AAV (n=411) and MPO-AAV (n=162) versus controls (n=1595). Disease-associated genetic variants were analysed for association with organ involvement, age at diagnosis and relapse, respectively. RESULTS: PR3-AAV was significantly associated with both HLA-DPB1*04:01 and rs1042335 at the HLA-DPB1 locus, also after stepwise conditional analysis. MPO-AAV was significantly associated with HLA-DRB1*04:04. Neither carriage of HLA-DPB1*04:01 alleles in PR3-AAV nor of HLA-DRB1*04:04 alleles in MPO-AAV were associated with organ involvement, age at diagnosis or relapse. CONCLUSIONS: The association to the HLA region was distinct in Scandinavian cases with MPO-AAV compared with cases of East Asian descent. In PR3-AAV, the two separate signals of association to the HLD-DPB1 region mediate potentially different functional effects.

Genome-wide association study identifies new loci associated with OCD.

To date, four genome-wide association studies (GWAS) of obsessive-compulsive disorder (OCD) have been published, reporting a high single-nucleotide polymorphism (SNP)-heritability of 28% but finding only one significant SNP. A substantial increase in sample size will likely lead to further identification of SNPs, genes, and biological pathways mediating the susceptibility to OCD. We conducted a GWAS meta-analysis with a 2-3-fold increase in case sample size (OCD cases: N = 37,015, controls: N = 948,616) compared to the last OCD GWAS, including six previously published cohorts (OCGAS, IOCDF-GC, IOCDF-GC-trio, NORDiC-nor, NORDiC-swe, and iPSYCH) and unpublished self-report data from 23andMe Inc. We explored the genetic architecture of OCD by conducting gene-based tests, tissue and celltype enrichment analyses, and estimating heritability and genetic correlations with 74 phenotypes. To examine a potential heterogeneity in our data, we conducted multivariable GWASs with MTAG. We found support for 15 independent genome-wide significant loci (14 new) and 79 protein-coding genes. Tissue enrichment analyses implicate multiple cortical regions, the amygdala, and hypothalamus, while cell type analyses yielded 12 cell types linked to OCD (all neurons). The SNP-based heritability of OCD was estimated to be 0.08. Using MTAG we found evidence for specific genetic underpinnings characteristic of different cohort-ascertainment and identified additional significant SNPs. OCD was genetically correlated with 40 disorders or traits-positively with all psychiatric disorders and negatively with BMI, age at first birth and multiple autoimmune diseases. The GWAS meta-analysis identified several biologically informative genes as important contributors to the aetiology of OCD. Overall, we have begun laying the groundwork through which the biology of OCD will be understood and described.

Genome wide association study in Swedish Labrador retrievers identifies genetic loci associated with hip dysplasia and body weight.

Genome wide association studies (GWAS) have been utilized to identify genetic risk loci associated with both simple and complex inherited disorders. Here, we performed a GWAS in Labrador retrievers to identify genetic loci associated with hip dysplasia and body weight. Hip dysplasia scores were available for 209 genotyped dogs. We identified a significantly associated locus for hip dysplasia on chromosome 24, with three equally associated SNPs (p = 4.3 × 10-7) in complete linkage disequilibrium located within NDRG3, a gene which in humans has been shown to be differentially expressed in osteoarthritic joint cartilage. Body weight, available for 85 female dogs, was used as phenotype for a second analysis. We identified two significantly associated loci on chromosome 10 (p = 4.5 × 10-7) and chromosome 31 (p = 2.5 × 10-6). The most associated SNPs within these loci were located within the introns of the PRKCE and CADM2 genes, respectively. PRKCE has been shown to play a role in regulation of adipogenesis whilst CADM2 has been associated with body weight in multiple human GWAS. In summary, we identified credible candidate loci explaining part of the genetic inheritance for hip dysplasia and body weight in Labrador retrievers with strong candidate genes in each locus previously implicated in the phenotypes investigated.

A survey of ficolin-3 activity in Systemic Lupus Erythematosus reveals a link to hematological disease manifestations and autoantibody profile.

The complement system plays a central role in the pathogenesis of Systemic Lupus Erythematosus (SLE), but most studies have focused on the classical pathway. Ficolin-3 is the main initiator of the lectin pathway of complement in humans, but its role in systemic autoimmune disease has not been conclusively determined. Here, we combined biochemical and genetic approaches to assess the contribution of ficolin-3 to SLE risk and disease manifestations. Ficolin-3 activity was measured by a functional assay in serum or plasma samples from Swedish SLE patients (n = 786) and controls matched for age and sex (n = 566). Genetic variants in an extended 300 kb genomic region spanning the FCN3 locus were analyzed for their association with ficolin-3 activity and SLE manifestations in a Swedish multicenter cohort (n = 985). Patients with ficolin-3 activity in the highest tertile showed a strong enrichment in an SLE cluster defined by anti-Sm/DNA/nucleosome antibodies (OR 3.0, p < 0.001) and had increased rates of hematological disease (OR 1.4, p = 0.078) and lymphopenia (OR = 1.6, p = 0.039). Genetic variants associated with low ficolin-3 activity mapped to an extended haplotype in high linkage disequilibrium upstream of the FCN3 gene. Patients carrying the lead genetic variant associated with low ficolin-3 activity had a lower frequency of hematological disease (OR 0.67, p = 0.018) and lymphopenia (OR 0.63, p = 0.031) and fewer autoantibodies (p = 0.0019). Loss-of-function variants in the FCN3 gene were not associated with SLE, but four (0.5 %) SLE patients developed acquired ficolin-3 deficiency where ficolin-3 activity in serum was depleted following diagnosis of SLE. Taken together, our results provide genetic and biochemical evidence that implicate the lectin pathway in hematological SLE manifestations. We also identify lectin pathway activation through ficolin-3 as a factor that contributes to the autoantibody response in SLE.

A Burden of Rare Copy Number Variants in Obsessive-Compulsive Disorder.

Current genetic research on obsessive-compulsive disorder (OCD) supports contributions to risk specifically from common single nucleotide variants (SNVs), along with rare coding SNVs and small insertion-deletions (indels). The contribution to OCD risk from large, rare copy number variants (CNVs), however, has not been formally assessed at a similar scale. Here we describe an analysis of rare CNVs called from genotype array data in 2,248 deeply phenotyped OCD cases and 3,608 unaffected controls from Sweden and Norway. We found that in general cases carry an elevated burden of large (>30kb, at least 15 probes) CNVs (OR=1.12, P=1.77×10-3). The excess rate of these CNVs in cases versus controls was around 0.07 (95% CI 0.02-0.11, P=2.58×10-3). This signal was largely driven by CNVs overlapping protein-coding regions (OR=1.19, P=3.08×10-4), particularly deletions impacting loss-of-function intolerant genes (pLI>0.995, OR=4.12, P=2.54×10-5). We did not identify any specific locus where CNV burden was associated with OCD case status at genome-wide significance, but we noted non-random recurrence of CNV deletions in cases (permutation P = 2.60×10-3). In cases where sufficient clinical data were available (n=1612) we found that carriers of neurodevelopmental duplications were more likely to have comorbid autism (P<0.001), and that carriers of deletions overlapping neurodevelopmental genes had lower treatment response (P=0.02). The results demonstrate a contribution of large, rare CNVs to OCD risk, and suggest that studies of rare coding variation in OCD would have increased power to identify risk genes if this class of variation were incorporated into formal tests.

Distinct HLA associations with autoantibody-defined subgroups in idiopathic inflammatory myopathies.

BACKGROUND: In patients with idiopathic inflammatory myopathies (IIM), autoantibodies are associated with specific clinical phenotypes suggesting a pathogenic role of adaptive immunity. We explored if autoantibody profiles are associated with specific HLA genetic variants and clinical manifestations in IIM. METHODS: We included 1348 IIM patients and determined the occurrence of 14 myositis-specific or -associated autoantibodies. We used unsupervised cluster analysis to identify autoantibody-defined subgroups and logistic regression to estimate associations with clinical manifestations, HLA-DRB1, HLA-DQA1, HLA-DQB1 alleles, and amino acids imputed from genetic information of HLA class II and I molecules. FINDINGS: We identified eight subgroups with the following dominant autoantibodies: anti-Ro52, -U1RNP, -PM/Scl, -Mi2, -Jo1, -Jo1/Ro52, -TIF1γ or negative for all analysed autoantibodies. Associations with HLA-DRB1∗11, HLA-DRB1∗15, HLA-DQA1∗03, and HLA-DQB1∗03 were present in the anti-U1RNP-dominated subgroup. HLA-DRB1∗03, HLA-DQA1∗05, and HLA-DQB1∗02 alleles were overrepresented in the anti-PM/Scl and anti-Jo1/Ro52-dominated subgroups. HLA-DRB1∗16, HLA-DRB1∗07 alleles were most frequent in anti-Mi2 and HLA-DRB1∗01 and HLA-DRB1∗07 alleles in the anti-TIF1γ subgroup. The HLA-DRB1∗13, HLA-DQA1∗01 and HLA-DQB1∗06 alleles were overrepresented in the negative subgroup. Significant signals from variations in class I molecules were detected in the subgroups dominated by anti-Mi2, anti-Jo1/Ro52, anti-TIF1γ, and the negative subgroup. INTERPRETATION: Distinct HLA class II and I associations were observed for almost all autoantibody-defined subgroups. The associations support autoantibody profiles use for classifying IIM which would likely reflect underlying pathogenic mechanisms better than classifications based on clinical symptoms and/or histopathological features. FUNDING: See a detailed list of funding bodies in the Acknowledgements section at the end of the manuscript.

Genome sequencing of 2000 canids by the Dog10K consortium advances the understanding of demography, genome function and architecture.

BACKGROUND: The international Dog10K project aims to sequence and analyze several thousand canine genomes. Incorporating 20 × data from 1987 individuals, including 1611 dogs (321 breeds), 309 village dogs, 63 wolves, and four coyotes, we identify genomic variation across the canid family, setting the stage for detailed studies of domestication, behavior, morphology, disease susceptibility, and genome architecture and function. RESULTS: We report the analysis of > 48 M single-nucleotide, indel, and structural variants spanning the autosomes, X chromosome, and mitochondria. We discover more than 75% of variation for 239 sampled breeds. Allele sharing analysis indicates that 94.9% of breeds form monophyletic clusters and 25 major clades. German Shepherd Dogs and related breeds show the highest allele sharing with independent breeds from multiple clades. On average, each breed dog differs from the UU_Cfam_GSD_1.0 reference at 26,960 deletions and 14,034 insertions greater than 50 bp, with wolves having 14% more variants. Discovered variants include retrogene insertions from 926 parent genes. To aid functional prioritization, single-nucleotide variants were annotated with SnpEff and Zoonomia phyloP constraint scores. Constrained positions were negatively correlated with allele frequency. Finally, the utility of the Dog10K data as an imputation reference panel is assessed, generating high-confidence calls across varied genotyping platform densities including for breeds not included in the Dog10K collection. CONCLUSIONS: We have developed a dense dataset of 1987 sequenced canids that reveals patterns of allele sharing, identifies likely functional variants, informs breed structure, and enables accurate imputation. Dog10K data are publicly available.

Using evolutionary constraint to define novel candidate driver genes in medulloblastoma.

Current knowledge of cancer genomics remains biased against noncoding mutations. To systematically search for regulatory noncoding mutations, we assessed mutations in conserved positions in the genome under the assumption that these are more likely to be functional than mutations in positions with low conservation. To this end, we use whole-genome sequencing data from the International Cancer Genome Consortium and combined it with evolutionary constraint inferred from 240 mammals, to identify genes enriched in noncoding constraint mutations (NCCMs), mutations likely to be regulatory in nature. We compare medulloblastoma (MB), which is malignant, to pilocytic astrocytoma (PA), a primarily benign tumor, and find highly different NCCM frequencies between the two, in agreement with the fact that malignant cancers tend to have more mutations. In PA, a high NCCM frequency only affects the BRAF locus, which is the most commonly mutated gene in PA. In contrast, in MB, >500 genes have high levels of NCCMs. Intriguingly, several loci with NCCMs in MB are associated with different ages of onset, such as the HOXB cluster in young MB patients. In adult patients, NCCMs occurred in, e.g., the WASF-2/AHDC1/FGR locus. One of these NCCMs led to increased expression of the SRC kinase FGR and augmented responsiveness of MB cells to dasatinib, a SRC kinase inhibitor. Our analysis thus points to different molecular pathways in different patient groups. These newly identified putative candidate driver mutations may aid in patient stratification in MB and could be valuable for future selection of personalized treatment options.

Relation between HLA and copy number variation of steroid 21-hydroxylase in a Swedish cohort of patients with autoimmune Addison's disease.

OBJECTIVE: Autoantibodies against the adrenal enzyme 21-hydroxylase is a hallmark manifestation in autoimmune Addison's disease (AAD). Steroid 21-hydroxylase is encoded by CYP21A2, which is located in the human leucocyte antigen (HLA) region together with the highly similar pseudogene CYP21A1P. A high level of copy number variation is seen for the 2 genes, and therefore, we asked whether genetic variation of the CYP21 genes is associated with AAD. DESIGN: Case-control study on patients with AAD and healthy controls. METHODS: Using next-generation DNA sequencing, we estimated the copy number of CYP21A2 and CYP21A1P, together with HLA alleles, in 479 Swedish patients with AAD and autoantibodies against 21-hydroxylase and in 1393 healthy controls. RESULTS: With 95% of individuals carrying 2 functional 21-hydroxylase genes, no difference in CYP21A2 copy number was found when comparing patients and controls. In contrast, we discovered a lower copy number of the pseudogene CYP21A1P among AAD patients (P = 5 × 10-44), together with associations of additional nucleotide variants, in the CYP21 region. However, the strongest association was found for HLA-DQB1*02:01 (P = 9 × 10-63), which, in combination with the DRB1*04:04-DQB1*03:02 haplotype, imposed the greatest risk of AAD. CONCLUSIONS: We identified strong associations between copy number variants in the CYP21 region and risk of AAD, although these associations most likely are due to linkage disequilibrium with disease-associated HLA class II alleles.

Mutations in the CYP27B1 gene cause vitamin D dependent rickets in pugs.

Rickets is a disorder of bone development and can be the result of either dietary or genetic causes. Here, related pugs from 2 litters were included. Three pugs had clinical signs including, lameness, bone deformities, and dyspnea. One other pug was found dead. Radiographs of 2 affected pugs, 5 and 6 months old, showed generalized widening, and irregular margination of the physes of both the appendicular and the axial skeleton with generalized decrease in bone opacity and bulbous swelling of the costochondral junctions. Two pugs had low serum calcium and 1,25 (OH)2 D3 concentrations. Test results further indicated secondary hyperparathyroidism with adequate concentrations of 25-hydroxyvitamin D. Necropsy revealed tongue-like projections of cartilage extending into the metaphysis consistent with rickets, loss of metaphyseal mineralization and lung pathology. Vitamin D-dependent rickets was diagnosed. A truncating mutation in the 1α-hydroxylase gene (CYP27B1) was identified by genome sequence analysis of the pugs with VDDR type 1A. Vitamin D-dependent rickets type 1A can occur in young pugs, and if left untreated is a life-threatening condition. Early medical intervention can reverse clinical signs and should be instituted as soon as possible.

Insights into mammalian TE diversity through the curation of 248 genome assemblies.

We examined transposable element (TE) content of 248 placental mammal genome assemblies, the largest de novo TE curation effort in eukaryotes to date. We found that although mammals resemble one another in total TE content and diversity, they show substantial differences with regard to recent TE accumulation. This includes multiple recent expansion and quiescence events across the mammalian tree. Young TEs, particularly long interspersed elements, drive increases in genome size, whereas DNA transposons are associated with smaller genomes. Mammals tend to accumulate only a few types of TEs at any given time, with one TE type dominating. We also found association between dietary habit and the presence of DNA transposon invasions. These detailed annotations will serve as a benchmark for future comparative TE analyses among placental mammals.

The contribution of historical processes to contemporary extinction risk in placental mammals.

Species persistence can be influenced by the amount, type, and distribution of diversity across the genome, suggesting a potential relationship between historical demography and resilience. In this study, we surveyed genetic variation across single genomes of 240 mammals that compose the Zoonomia alignment to evaluate how historical effective population size (Ne) affects heterozygosity and deleterious genetic load and how these factors may contribute to extinction risk. We find that species with smaller historical Ne carry a proportionally larger burden of deleterious alleles owing to long-term accumulation and fixation of genetic load and have a higher risk of extinction. This suggests that historical demography can inform contemporary resilience. Models that included genomic data were predictive of species' conservation status, suggesting that, in the absence of adequate census or ecological data, genomic information may provide an initial risk assessment.

Mammalian evolution of human cis-regulatory elements and transcription factor binding sites.

Understanding the regulatory landscape of the human genome is a long-standing objective of modern biology. Using the reference-free alignment across 241 mammalian genomes produced by the Zoonomia Consortium, we charted evolutionary trajectories for 0.92 million human candidate cis-regulatory elements (cCREs) and 15.6 million human transcription factor binding sites (TFBSs). We identified 439,461 cCREs and 2,024,062 TFBSs under evolutionary constraint. Genes near constrained elements perform fundamental cellular processes, whereas genes near primate-specific elements are involved in environmental interaction, including odor perception and immune response. About 20% of TFBSs are transposable element-derived and exhibit intricate patterns of gains and losses during primate evolution whereas sequence variants associated with complex traits are enriched in constrained TFBSs. Our annotations illuminate the regulatory functions of the human genome.

A genomic timescale for placental mammal evolution.

The precise pattern and timing of speciation events that gave rise to all living placental mammals remain controversial. We provide a comprehensive phylogenetic analysis of genetic variation across an alignment of 241 placental mammal genome assemblies, addressing prior concerns regarding limited genomic sampling across species. We compared neutral genome-wide phylogenomic signals using concatenation and coalescent-based approaches, interrogated phylogenetic variation across chromosomes, and analyzed extensive catalogs of structural variants. Interordinal relationships exhibit relatively low rates of phylogenomic conflict across diverse datasets and analytical methods. Conversely, X-chromosome versus autosome conflicts characterize multiple independent clades that radiated during the Cenozoic. Genomic time trees reveal an accumulation of cladogenic events before and immediately after the Cretaceous-Paleogene (K-Pg) boundary, implying important roles for Cretaceous continental vicariance and the K-Pg extinction in the placental radiation.

​Comparative genomics of Balto, a famous historic dog, captures lost diversity of 1920s sled dogs.

We reconstruct the phenotype of Balto, the heroic sled dog renowned for transporting diphtheria antitoxin to Nome, Alaska, in 1925, using evolutionary constraint estimates from the Zoonomia alignment of 240 mammals and 682 genomes from dogs and wolves of the 21st century. Balto shares just part of his diverse ancestry with the eponymous Siberian husky breed. Balto's genotype predicts a combination of coat features atypical for modern sled dog breeds, and a slightly smaller stature. He had enhanced starch digestion compared with Greenland sled dogs and a compendium of derived homozygous coding variants at constrained positions in genes connected to bone and skin development. We propose that Balto's population of origin, which was less inbred and genetically healthier than that of modern breeds, was adapted to the extreme environment of 1920s Alaska.

The functional and evolutionary impacts of human-specific deletions in conserved elements.

Conserved genomic sequences disrupted in humans may underlie uniquely human phenotypic traits. We identified and characterized 10,032 human-specific conserved deletions (hCONDELs). These short (average 2.56 base pairs) deletions are enriched for human brain functions across genetic, epigenomic, and transcriptomic datasets. Using massively parallel reporter assays in six cell types, we discovered 800 hCONDELs conferring significant differences in regulatory activity, half of which enhance rather than disrupt regulatory function. We highlight several hCONDELs with putative human-specific effects on brain development, including HDAC5, CPEB4, and PPP2CA. Reverting an hCONDEL to the ancestral sequence alters the expression of LOXL2 and developmental genes involved in myelination and synaptic function. Our data provide a rich resource to investigate the evolutionary mechanisms driving new traits in humans and other species.

Integrating gene annotation with orthology inference at scale.

Annotating coding genes and inferring orthologs are two classical challenges in genomics and evolutionary biology that have traditionally been approached separately, limiting scalability. We present TOGA (Tool to infer Orthologs from Genome Alignments), a method that integrates structural gene annotation and orthology inference. TOGA implements a different paradigm to infer orthologous loci, improves ortholog detection and annotation of conserved genes compared with state-of-the-art methods, and handles even highly fragmented assemblies. TOGA scales to hundreds of genomes, which we demonstrate by applying it to 488 placental mammal and 501 bird assemblies, creating the largest comparative gene resources so far. Additionally, TOGA detects gene losses, enables selection screens, and automatically provides a superior measure of mammalian genome quality. TOGA is a powerful and scalable method to annotate and compare genes in the genomic era.

Three-dimensional genome rewiring in loci with human accelerated regions.

Human accelerated regions (HARs) are conserved genomic loci that evolved at an accelerated rate in the human lineage and may underlie human-specific traits. We generated HARs and chimpanzee accelerated regions with an automated pipeline and an alignment of 241 mammalian genomes. Combining deep learning with chromatin capture experiments in human and chimpanzee neural progenitor cells, we discovered a significant enrichment of HARs in topologically associating domains containing human-specific genomic variants that change three-dimensional (3D) genome organization. Differential gene expression between humans and chimpanzees at these loci suggests rewiring of regulatory interactions between HARs and neurodevelopmental genes. Thus, comparative genomics together with models of 3D genome folding revealed enhancer hijacking as an explanation for the rapid evolution of HARs.

Leveraging base-pair mammalian constraint to understand genetic variation and human disease.

Thousands of genomic regions have been associated with heritable human diseases, but attempts to elucidate biological mechanisms are impeded by an inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function, agnostic to cell type or disease mechanism. Single-base phyloP scores from 240 mammals identified 3.3% of the human genome as significantly constrained and likely functional. We compared phyloP scores to genome annotation, association studies, copy-number variation, clinical genetics findings, and cancer data. Constrained positions are enriched for variants that explain common disease heritability more than other functional annotations. Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.